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PREFACE 


This  report  was  prepared  by  Paul  W.  Richmond,  Mechanical  Engineer,  George  L.  Blaisdell, 
Research  Civil  Engineer,  both  of  the  Applied  Research  Branch,  Experimental  Engineering  Division, 
U.S.  Army  Cold  Regions  Research  and  Engineering  Laboratory,  and  Charles  E.  Green,  Research  Civil 
Engineer,  Mobility  Investigations  Group,  Mobility  Systems  Division,  U.S.  Army  Waterways  Experiment 
Station  (WES).  The  report  documents  some  of  the  efforts  expended  in  a  joint  study  conducted  by 
CRREL  and  WES.  The  Keweenaw  Research  Center  (KRC)  of  Michigan  Technological  University 
was  contracted  with  to  provide  services  and  expertise  for  this  project  as  well.  Funding  for  this  report 
was  provided  by  DA  Project  4A762784AT42,  Cold  Regions  Engineering  Technology,  Work  Unit  CS/ 
040,  Wheels  vs  Tracks  in  Winter. 

Technical  review  of  this  report  was  graciously  provided  by  Sally  A.  Shoop  (CRREL)  and  Russell 
Alger  (KRC). 

The  authors  express  their  appreciation  to  the  personnel  at  CRREL,  WES  and  KRC  who  worked 
on  this  project.  Particular  thanks  go  to  Stephen  Decato  and  Andrew  Sunderlund  of  CRREL  and  to 
Russell  Alger  of  KRC. 

The  contents  of  this  report  are  not  to  be  used  for  advertising  or  promotional  purposes.  Citation  of 
brand  names  does  not  constitute  an  official  endorsement  or  approval  of  the  use  of  such  commercial 
products. 
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length  of  tire  or  track  in  contact  with  undefonned  snow 
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average  width  of  the  deformed  snow  under  a  wheel  or  track 
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draw  bar  pull 
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tire  or  track  contact  pressure 
tire  radius 

hard  surface  motion  resistance 

external  resistance  attributable  to  snow  compaction 

total  motion  resistance 

gross  traction 

towed  motion  resistance 

net  traction 

vehicle  weight 

vehicle  sinkage 

maximum  sinkage  for  a  vehicle 
theoretical  final  density 
initial  density 


Abbreviations 

CTV  CRREL  Instrumented  Vehicle 
HEMTT  Heavy  Expanded  Mobility  Tactical  Truck 
HMMWV  High-Mobility  Multipurpose  Wheeled  Vehicle 
KRC  Keweenaw  Research  Center 
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Wheels  and  TVacks  in  Snow 

Second  Validation  Study  of  the  CRREL  Shallow  Snow  Mobility  Model 

PAUL  W.  RICHMOND,  GEORGE  L.  BLAISDELL 
AND  CHARLES  E.  GREEN 


INTRODUCTION 

During  the  winters  of  1 988  and  1 989,  a  winter  mobil¬ 
ity  study  was  jointly  conducted  by  WES  and  CRREL. 
These  studies,  part  of  the  larger  U.  S.  Army  Wheels / 
Tracks  program,  were  to  be  used  to  compare  predictions 
of  the  CRREL  shallow  snow  mobility  model  with  actual 
snow  mobility  data  for  a  wide  variety  of  vehicles.  Some 
of  the  1988  results  are  reported  in  Blaisdell  et  al.  (1990) 
and  Green  and  Blaisdell  (in  press).  The  winter  phase  of 
the  WheelsITracks  study  encompassed  two  winter  field 
seasons;  this  report  presents  the  results  of  the  second 
winter  field  season  (1989). 

The  major  accomplishment  of  the  first  (1988)  field 
study  was  the  development  of  a  new  traction  algorithm, 
which  was  incorporated  into  the  second  version  of  the 
shallow  snow  mobility  model  (SSM2.0).  This  traction 
equation  was  based  on  shear  failure  of  snow  via  the  Mohr- 
Coulomb  criterion  and  used  test  data  to  arrive  at  a  best- 
fit  failure  curve.  The  curve  was  found  to  predict  traction 
well  for  all  snow  densities  because  tractive  effort  in  most 
cases  occurs  on  similar  snow  conditions  (compacted 
snow  with  a  density  of  approximately  550  kg/m3)  re¬ 
gardless  of  the  initial  snow  conditions.  Additionally,  the 
equation  was  developed  using  data  from  wheeled  vehicles 
equipped  with  state-of-the-art  tires  and  several  tracked 
vehicles.  Vehicle  motion  >  distance,  however,  was  not 
predicted  well  by  the  shallow  snow  model. 

Based  on  the  results  and  analysis  of  the  1 98 8  data,  we 
decided  that  the  final  winter  field  season  in  the  Wheels / 
Tracks  program  would  primarily  address  vehicle  motion 
resistance.  Limited  attention  would  be  given  to  the 
traction  aspect  of  mobility  and  this  would  be  oriented 
towards  removing  some  of  the  known  caveats  in  the 
SSM2.0  traction  algorithm. 

During  the  second  field  season,  tests  were  conducted 
primarily  during  January  1 989  with  limited  testing  con¬ 
tinuing  in  March  and  April  1989. 


BACKGROUND 

The  primary  goal  of  this  study  was  to  investigate 
vehicle  motion  resistance  in  snow  to  continue  validation 
of  the  CRREL  shallow  snow  mobility  model.  The  model 
is  based  on  theoretical  relationships  and  empirical  ex¬ 
pressions  developed  in  the  past  from  a  large,  but  scat¬ 
tered,  data  base. 

Briefly,  SSM2.0  uses  the  following  expressions  to 
calculate  shallow  snow  mobility: 

Net  traction:  Tn  =  Tg-Rs  ( 1 ) 

where  T&  is  gross  traction  and  Rs  is  the  external  motion 
resistance  attributable  to  snow  compaction.  Gross  trac¬ 
tion  (in  kilopascals)  can  be  estimated  fora  wide  range  of 
vehicles  using 

Tg  =  0.851  A/0-823  (2) 

where  N  is  the  normal  stress  under  a  tire  or  track  (in 
kilopascals).  Rs  is  primarily  determined  by  the  amount  a 
vehicle  sinks  in  the  snow;  SSM2.0  uses  the  following 
equations  to  estimate  sinkage  (z) 

Maximum  sinkage:  zmax  =  A  (1  —  (3) 

Pf 

where  h  is  the  depth  of  undisturbed  snow  (in  meters),  pf 
is  the  theoretical  final  density  in  the  rut  following  vehicle 
passage  (kilograms  per  cubic  meter,  see  Table  1  for 
values)  and  pf  is  the  initial  undisturbed  snow  density 
(kilograms  per  cubic  meter). 

Both  SSM1.0  and  SSM2.0  were  described  fully  by 
Blaisdell  et  al.  ( 1 990);  we  repeat  a  detailed  description  of 
SSM2.0  in  Appendix  A  for  reference,  and  refer  the 
reader  to  Appendix  A  for  the  equation  describing  Rr 


Table  1.  Final  snow  densities  used 
in  SSM2.0. 


Maximum  ground  pressure  Final  density 

(kPa)  (kg/in3) 


<210 

500 

211-350 

550 

351-700 

600 

>701 

650 

FIELD  EXPERIMENTS 

The  field  experiments  described  in  this  report  were 
carried  out  at  Keweenaw  Research  Center  (KRC)  lo¬ 
cated  at  the  Houghton  County  Airpark,  Michigan  (KRC 
is  located  on  the  Keweenaw  Peninsula  of  Lake  Supe¬ 
rior).  Tests  were  conducted  by  personnel  from  WES, 
KRC  and  CRREL.  Two  types  of  mobility  tests  were 
conducted — traction  tests  in  various  snow  conditions 
and  resistance  tests  in  undisturbed  snow. 

CRREL  and  KRC  personnel  conducted  mobility  tests 
using  the  CRREL  Instrumented  Vehicle  (CIV).  The 
CIV,  which  is  fully  described  by  Blaisdell  (1983),  is 
based  on  a  1 977  Jeep  Cherokee  with  modifications  to  its 
braking  and  driving  components  to  accommodate  typi¬ 
cal  mobility  tests.  In  addition  to  its  onboard  computer- 
based  data  acquisition  system,  the  vehicle  contains  a 


number  of  transducers  for  force  and  speed  measure¬ 
ments.  In  previous  tests  the  CIV  has  produced  results 
(traction  data)  that  agreed  very  well  with  data  obtained 
from  larger  vehicles  with  similar  ground  contact  pres¬ 
sures  (Blaisdell  et  al .,  1 990)  in  side  by  side  tests.  For  this 
reason  the  CIV  was  used  again  this  field  season  to  extend 
the  winter  mobility  data  base.  The  CIV  was  equipped 
with  several  different  tires  during  this  study  to  examine 
the  effects  of  tire  parameters  on  traction  and  resistance. 

The  standard  or  control  tire  was  a  Michelin  XCH4  all- 
season  mud  and  snow  tire  and  was  tested  at  inflation 
pressures  of  179  and  103  kPa  (26  and  15  lb/in.2).  During 
the  traction  tests  two  additional  tires  were  used:  an 
NDCC  (Non-Directional  Cross  Country)  tire,  which  is 
the  old  standard  military  tire,  and  another  Michelin 
XCH4  that  had  its  tread  buffed  off.  The  NDCC  tire  was 
to  represent  a  tire  that  wasn’t  up  to  the  latest  standards  in 
tire  design,  and  the  buffed  Michelin  was  used  to  deter¬ 
mine  the  effect  of  tire  tread,  as  compared  with  the 
unmodified  control  tire. 

For  the  resistance  tests,  the  control  tire  and  two  other 
tires  were  used.  These  two  additional  tires  were  not  used 
for  traction  tests  and  were  chosen  for  resistance  testing 
based  solely  on  their  width  and  availability.  The  Firestone 
T 1 45/80  is  a  temporary  spare  tire  with  a  maximum  width 
of  0. 1 56  m  (6. 1  in.)  and  the  Goodyear  Eagle  P225/60R 1 5 
is  a  low  profile  “street-rod”  tire  with  a  maximum  width 
of  0.274  m  (10.8  in.).  The  tires  used  with  the  CIV  and 
their  characteristics  are  given  in  Table  2. 


Table  2.  CIV  tire  data. 


Inflation  Contact  Contact 

Tire  pressure  area*  Radius *  Width**  length  Tire 


nomenclature 

(kPa) 

(n?) 

(m) 

< m ) 

(m) 

code 

Michelin 

179 

0.0412 

0.375 

0.260 

0.2512 

A 

LT235XCH4 

103 

0.0635 

0.375 

0.279 

0.3825 

A 

Michelin  (buffed)™ 

179 

0.0443 

0.37 

0.260 

0.253 

a 

LT235XCH4 

103 

0.0586 

0.37 

0.272 

0.335 

a 

NDCC 

234 

0.0238 

0.39 

0.215 

0.1993 

B 

700-16LW 

138 

0.029 

0.39 

0.217 

0.2377 

B 

Firestone 

414 

0.0204 

0.325 

0.156 

0.187 

C 

T145/80 

Goodyear  Eagle 

248 

0.0319 

0.35 

0.274 

0.1636 

D 

P225/60R15 

Goodyear  Tiempo 

179 

0.028 

0.356 

0.254 

0.2073 

E 

P225/75R15 

103 

0.034 

0.349 

0.267 

0.2564 

*  Hard  surface  contact  area. 

*  Radius  of  undeformed  part  of  lire. 

**  Maximum  deformed  width. 

ft  The  tread  was  buffed  off  to  below  the  wear  bars. 
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Table  3.  Tire  and  track  data  for  selected  vehicles. 


Tire 

Inflation 

pressure 

Contact 

area* 

Radius * 

Width** 

Contact 

length 

Vehicle 

nomenclature 

(kPa) 

(n?) 

<m) 

(mi 

(m) 

HMMWV 

37.00  x  12.5  R16.5 

138/152 

0.074 

0.429 

0.33 

0.247 

36.00  x  12.5  LT 

138 

0.425 

HEMTT 

16.0  R20 

241/276 

0.149 

0.617 

0.475 

0.429 

139/207 

0.171 

0.589 

0.483 

0.472 

LAV25** 

12.50/75  R  20XL 

207 

0.100 

0.445 

0.378 

0.368 

103 

0.141 

0.414 

0.343 

0.518 

1 1.00  R  16XL 

290 

0.580 

0.434 

0.314 

_ 

165 

0.102 

0.417 

0.332 

— 

SUSV 

Track 

— 

1.18 

32.5° 

0.6096 

3.7592 

Bradley** 

Track 

— 

2.09 

25° 

0.533 

3.920 

M113A1** 

Track 

— 

1.02 

21° 

0.381 

2.667 

M601A1 

Track 

— 

2.8 

35° 

0.66 

4.248 

*  Per  tire  on  track. 

*  Radius  for  wheeled  vehicles,  entrance  angle  (degrees)  for  tracked  vehicles. 
**  Maximum  deformed  width. 

**  From  1988  tests. 


Personnel  from  WES  conducted  mobility  tests  using 
the  following  group  of  military  vehicles: 

1.  M988 — High-Mobility  Multipurpose  Wheeled 
Vehicle  (HMMWV),  4x4,  equipped  with  Michelin 
37xl2.5R16.5LT  tires. 

2.  M977 — Heavy  Expanded  Mobility  Tactical  truck 
(HEMTT),  8x8,  equipped  with  Michelin  16.0R20  tires. 

3.  M973— Small  Unit  Support  Vehicle  (SUSV),  ar¬ 
ticulated,  tracked. 

4.  M60A1 — Main  battle  tank,  tracked. 

These  vehicles  were  chosen  to  span  the  full  range  of 
typical  ground  vehicle  contact  pressures  and  to  represent 
a  cross  section  of  current  military  vehicles.  The  charac¬ 
teristics  of  these  vehicles,  as  well  as  those  used  during 
the  1988  field  season,  are  given  in  Table  3. 

Test  procedures 

The  test  procedures  followed  typical  mobility  field 
studies,  in  that  measurements  of  net  traction  Tn  and  total 
motion  resistance  /?,  were  made  with  each  vehicle  under 
varying  snow  conditions.  Although  we  wanted  to  con¬ 
duct  tests  using  all  the  vehicles  on  each  day,  the  shortness 
of  the  available  field  time  and  the  lack  of  appropriate 
snow  falls  precluded  this.  All  of  the  tests  were  done  in 
areas  that  had  a  packed  snow  base,  with  the  exception  of 
some  traction  and  resistance  tests,  which  were  done  on 
a  snow-covered  area  that  was  underlaid  by  ice  (the  KRC 
ice  rink). 


CIV 

The  CIV’s  resistance  to  motion  is  measured  with  its 
rear  tires  driving  and  its  front  wheels  rolling  free.  Since 
the  triaxial  load  cells  are  located  just  inside  the  front 
wheels,  this  test  measures  the  total  amount  of  resistance 
felt  at  the  front  tires  only.  Motion  resistance  is  first 
established  on  a  level,  undeformable  surface.  Measure¬ 
ments  of  hard  surface  resistance  Rh  are  obtained  for  each 
tire  type  and  selected  inflation  pressure.  By  convention, 
motion  resistance  tests  are  conducted  at  a  vehicle  speed 
of  8  km/hr  (5  mi/hr),  and  it  is  known  that  resistance 
values  are  independent  of  moderate  variations  (±  3.2 
km/hr  [±  2  mi/hr])  in  vehicle  speed.  Variations  in  Rh 
values  between  tire  types  are  the  result  of  differences  in 
the  forces  necessary  for  tire  flexing  and  can  be  attributed 
to  differences  in  their  design. 

The  external  vehicle  motion  resistance  in  a  snow 
cover  (Rs)  is  calculated  by  measuring  the  motion  resis¬ 
tance  (/?,)  in  the  test  area  using  the  above  procedure  and 
subtracting  the  hard  surface  motion  resistance  (/?h). 

Traction  is  measured  by  accelerating  the  front 
(driving)  wheels  while  braking  the  rear  wheels  to  hold 
the  vehicle  speed  constant  at  8  km/hr  ±  1  km/hr  (5  ±  1/2 
mi/hr).  In  this  manner  the  front  tires  are  driven  through 
a  wide  range  of  slip  values,  starting  at  zero  slip.  A  plot  of 
measured  Tn  vs  differential  interface  velocity  (DIV)  is 
used  to  obtain  tractive  effort.  The  Tn  value  is  taken  as  the 
average  tractive  force  in  a  window  centered  around  the 
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maximum  tractive  force  reading.  The  window  is  chosen 
to  represent  a  range  of  slip  values  that  could  reasonably 
be  maintained  by  a  vehicle  operator.  Gross  traction  Tg  is 
then  calculated  from  eq  1  for  each  pair  of  resistance  and 
traction  tests. 

Military  vehicles 

The  tractive  performance  of  the  military  vehicles  was 
determined  by  measuring  drawbar  pull  ( DBP )  and  towed 
motion  resistance  ( TMR ).  These  measurements  are  not 
exactly  the  same  as  the  Tn  and  /?,  obtained  with  the  CIV, 
since  these  were  standard  military  vehicles  and  not 
modified  for  research.  Drawbar  pull  tests  were  con¬ 
ducted  by  measuring  the  force  that  a  vehicle  can  exert  on 
a  cable  that  is  being  used  to  resist  vehicle  motion.  Thus 

Tn  =QBP  (4) 

n' 

where  /?’  is  the  number  of  driven  tires  or  tracks.  Tn 
represents  the  average  tractive  force  per  driven  wheel  or 
track. 

Motion  resistance  is  determined  by  measuring  the 
vehicle’s  resistance  to  towing.  Here, 

/?,  =  TMR.  (5) 

n 

where  n  is  the  number  of  tires  or  tracks  and  R{  is  the  total 
motion  resistance  on  an  average  tire  or  track. 

Our  procedure  for  measuring  DBP  was  as  follows.  A 
load  vehicle  of  approximately  the  same  size  and  perfor¬ 
mance  as  the  test  vehicle  was  selected  to  apply  a  resis¬ 
tance  to  the  test  vehicle.  A  steel  cable  0.016  m  in 
diameter  (0.6  in.)  and  15.3  m  in  length  (50  ft)  was 
connected  from  the  front  of  the  load  vehicle  to  a  load  cell 
attached  to  the  rear  pintle  hook  of  the  test  vehicle.  A 
string  payout  system  (fishing  reel)  for  measurement  of 
true  ground  distance  was  also  mounted  on  the  test 
vehicle.  Tachometers  were  mounted  on  the  drive  wheels 
or  sprockets  of  the  test  vehicle  and  used  to  measure 
wheel  or  track  travel  during  a  test. 

During  each  test,  the  test  vehicles  were  operated  in 
their  lowest  gear  and  at  optimum  engine  rpm,  yielding 
vehicle  speeds  between  3.2  and  8  km/hr  (2  and  5  mi/hr). 
The  vehicle  proceeded  into  the  test  lane  with  the  load 
vehicle  following,  the  cable  between  the  two  vehicles 
being  slack  and  unloaded.  The  load  vehicle  driver  gradu¬ 
ally  applied  load  to  the  test  vehicle  by  braking.  The  test 
sequence  proceeded  from  the  test  vehicle  initially  expe¬ 
riencing  a  no-load,  no-slip  condition  and  increased  up  to 
a  high-load,  high-slip  condition  or  a  power  limited 
condition  in  which  the  test  vehicle  could  not  maintain  the 
desired  track  or  tire  speed.  Forward  speed  and  wheel  or 
track  speed  were  maintained  (very  nearly)  constant  for 
several  seconds  of  steady-state  pull  measurements.  Ve¬ 


hicle  slip  was  calculated  from  a  magnetic  tape  record  by 
using  both  vehicle  travel  distance  and  wheel  travel 
distance.  The  vehicle  slip,  in  percent,  is  equal  to 

x  100  (6) 

dv 

where  c/t  =  apparent  distance  traveled  by  the  wheel  or 

track 

dv  =  distance  actually  traveled  by  the  vehicle. 

Continuous  measurements  were  made  in  this  manner 
until  a  sufficient  number  of  load  and  slip  combinations 
were  recorded  to  develop  a  curve  of  drawbar  pull-slip 
data  (usually  two  good  test  sequences  in  the  same  area). 
As  with  Tn  in  the  CIV  traction  test,  DBP  is  a  function  of 
slip  or  DIV.  Generally,  in  low  density  snow  (less  than 
500  kg/m3)  maximum  DBP  occurs  at  rates  of  slip  greater 
than  20%.  However,  efficiency  of  operation  is  inversely 
proportional  to  slip;  little  useful  work  is  being  done  at 
very  high  slip  rates.  A  slipof  20%  is  generally  used  as  the 
cut  off  for  power  efficiency.  Thus,  the  maximum  DBP 
that  occurred  in  the  vicinity  of  20%  slip  was  used  for  the 
calculation  of  Tn  (eq  4).  Equation  1 ,  again,  is  engaged  to 
determine  gross  traction.  This  procedure  agrees  with  the 
CIV  data  analysis  process,  which  averaged  the  upper 
15%  of  the  gross  traction  data. 

The  procedure  used  for  obtaining  the  TMR  of  the  test 
vehicle  in  each  test  area  was  to  tow  the  vehicle  (with  its 
transmission  in  neutral)  at  a  speed  of  approximately  3.2 
km/hr  (2  mi/hr).  After  each  traction  test,  the  vehicle  was 
steered  into  an  undisturbed  area  adjacent  to  the  traction 
test  area,  usually  in  a  position  straddling  the  ruts  of  the 
associated  traction  test.  The  load  vehicle  then  towed  the 
test  vehicle  forward  to  determine  the  TMR.  The  test 
proceeded  for  a  sufficient  distance  to  permit  the  load  cell 
readings  to  stabilize  and  to  be  recorded  on  magnetic  tape. 
An  average  value  during  the  stable  portion  of  the  record 
was  taken  as  TMR,  and  /?,  was  then  calculated  from  eq  5. 
Finally,  external  resistance  from  snow  compaction  Rs 
was  found  by  subtracting  /?h  from  /?,. 

Measurements  of  snow  characteristics  (depth,  den¬ 
sity,  temperature  and  sinkage)  were  obtained  at  each  test 
area  while  tests  were  being  conducted.  The  undisturbed 
densities  ranged  from  70  to  320  kg/m3  (4.4  to  20  lb/ft3), 
and  depths  ranged  from  3  to  30  cm  (1  to  12  in.).  Gener¬ 
ally,  each  test  condition  represented  snow  from  one  storm 
and  the  air  temperatures  were  well  below  freezing. 


RESULTS  AND  ANALYSIS 
Traction 

The  objective  of  traction  tests  during  this  series  of 
tests  was  to  obtain  traction  data  for  snow  and  tire  condi- 
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Table  4.  CIV  traction  data. 


Tire 

Normal  load 

Normal  stress 

Resistance 

Net  traction 

Gross  traction 

Average 

Date 

Code 

Pressure 

Left 

Right 

Left 

Right 

Left 

Right 

Left 

Right 

Left 

Right 

Traction 

Load 

(1989) 

(kPa) 

(N) 

(N) 

(VP a) 

(kPa) 

(N) 

(N) 

(N) 

(N) 

(kPa) 

(kPa) 

(kPa) 

(kPa) 

17  Jan* 

A 

179 

6876 

167 

128 

172 

2323 

2124 

59 

56 

57.6 

173.3 

A 

179 

7369 

6913 

179 

168 

128 

172 

1859 

1696 

48 

45 

46.8 

173.3 

A 

103 

7514 

6912 

118 

109 

338 

348 

2076 

1978 

38 

37 

37.3 

113.6 

A 

103 

7509 

6905 

118 

109 

338 

348 

1949 

2010 

36 

37 

36.6 

113.5 

a 

179 

7488 

6598 

169 

149 

128 

147 

2679 

2483 

63 

59 

61.4 

159.0 

a 

179 

7455 

6486 

168 

146 

128 

147 

2854 

2705 

67 

64 

65.8 

157.4 

a 

103 

7336 

6869 

125 

117 

276 

253 

2977 

2788 

56 

52 

53.7 

121.2 

a 

103 

7543 

6747 

129 

115 

276 

253 

2502 

2297 

47 

44 

45.5 

121.9 

B 

234 

7631 

7015 

321 

295 

79 

59 

1441 

1190 

64 

52 

58.2 

307.7 

B 

234 

7561 

6886 

318 

289 

79 

59 

1474 

1337 

65 

59 

62.0 

303.5 

B 

138 

7554 

6990 

238 

220 

179 

134 

1486 

1366 

52 

47 

49.8 

228.7 

B 

138 

7506 

7029 

236 

221 

179 

134 

1489 

1566 

52 

53 

52.9 

228.5 

1  Marf 

A 

179 

7312 

6725 

177 

163 

786 

713 

1645 

1421 

59 

52 

55.4 

170.4 

A 

179 

7184 

6559 

174 

159 

786 

713 

1474 

1470 

55 

53 

53.9 

166.8 

A 

103 

7251 

6672 

114 

105 

795 

818 

1693 

1491 

39 

36 

37.8 

109.6 

A 

103 

7368 

6647 

118 

106 

795 

818 

1662 

1465 

39 

37 

37.9 

112.1 

a 

179 

7483 

6886 

169 

155 

885 

863 

893 

923 

40 

40 

40.2 

162.2 

a 

179 

7293 

6617 

165 

149 

885 

863 

1466 

1107 

53 

44 

48.8 

157.0 

a 

103 

7171 

6618 

122 

113 

1055 

876 

1934 

1770 

51 

45 

48.1 

117.7 

a 

103 

7188 

6419 

123 

110 

1055 

876 

1832 

1697 

49 

44 

46.6 

116.1 

B 

234 

7404 

6748 

311 

284 

578 

624 

1697 

1173 

96 

75 

85.6 

297.3 

B 

234 

7332 

6708 

308 

282 

578 

624 

1762 

1544 

98 

91 

94.7 

295.0 

B 

138 

7332 

6719 

231 

211 

657 

696 

1616 

1241 

71 

61 

66.2 

220.9 

B 

138 

7295 

6710 

229 

211 

657 

696 

1512 

1147 

68 

58 

63.1 

220.2 

6  Apr** 

A 

179 

6477 

6179 

157 

150 

842 

756 

1978 

1901 

68 

64 

66.2 

153.6 

A 

179 

6383 

6257 

155 

152 

842 

756 

2202 

2256 

74 

73 

73.3 

153.4 

A 

179 

6278 

6185 

152 

150 

842 

756 

2222 

2219 

74 

72 

73.1 

151.3 

A 

179 

6544 

6212 

159 

151 

149 

163 

2226 

2249 

57 

58 

57.8 

154.8 

A 

179 

6549 

6247 

159 

152 

149 

163 

2273 

2200 

58 

57 

57.8 

155.3 

*  On  17  January  the  tests  were  done  on  a  packed  and  groomed  snow  road,  the  snow  was  approximately  6  cm  deep  and  had  an  average  density  of  560  kg/m3. 
f  On  1  March  the  tests  were  done  on  undisturbed  dry  snow,  the  initial  snow  density  was  150  kg/m3  and  the  depth  ranged  from  10  to  20  cm. 

**On  6  April  the  snow  was  wet,  the  average  density  was  510  kg/m3  and  the  depth  was  2-7  cm,  the  last  two  tests  on  this  date  were  done  on  packed  wet  snow. 


tions  distinctly  different  from  those  given  in  Blaisdell  et 
al.  (1990),  which  were  used  to  develop  eq  2.  These  new 
conditions  were  to  be  used  to  extend  the  usefulness  of  eq 
2  from  undisturbed  snow  with  state-of-the-art  tires  to 
other  tires  and  snow  conditions;  additionally,  a  feel  for 
the  general  applicability  (robustness)  of  the  equation 
would  be  obtained.  Using  the  CIV,  we  examined  three 
additional  initial  snow  conditions  fortraction — wet  snow, 
wet  packed  snow  and  a  groomed  packed  snow  road.  The 
military  vehicles  were  tested  for  traction  on  the  packed 
snow  road  and  in  undisturbed  snow  overlaying  ice. 

CIV  traction 

The  traction  data  obtained  using  the  CIV,  sequen¬ 
tially  equipped  with  three  different  tire  types,  and  snow 
conditions  are  given  in  Table  4.  The  load  values  are  those 
measured  with  the  vehicle’s  load  cells.  The  stress  is 
obtained  by  dividing  the  load  by  the  hard  surface  contact 
areas  given  in  Table  2  for  each  tire  and  inflation  pressure. 
These  stress  values  are  plotted  in  Figure  1 ;  eq  2  (undis¬ 


turbed  snow,  underlain  by  packed  snow  or  frozen  soil)  is 
drawn  on  the  plots  for  reference.  There  is  good  agree¬ 
ment  between  the  traction  equation  and  the  data  under  all 
but  two  conditions.  The  performance  of  the  NDCC  tire 
on  packed  snow  is  over-predicted  by  eq  2,  and  the 
Michelin  XCH4’s  performance  is  under-predicted  on 
wet  snow.  The  performance  of  the  NDCC  tire  may  be 
explained  by  considering  tread  geometry  (deep,  widely 
spaced  lugs,  see  Appendix  B)  and  the  older  style  tire 
compound  and  carcass  design.  In  undisturbed  snow,  the 
lugs  of  the  NDCC  tire  were  able  to  dig  into  the  snow, 
developing  the  same  magnitude  of  interfacial  strength 
(traction)  as  the  state-of-the-art  tires.  However,  on  packed 
snow  the  less-than-optimal  tread  design  and  tire  com¬ 
pound  could  not  engage  the  same  amount  of  interfacial 
strength  as  observed  with  newer  type  of  tires.  In  wet 
snow  the  Michelin  XCH4  was  able  to  develop  more 
traction  than  expected.  This  may  be  ascribable  to  some 
uniqueness  of  the  tire’s  compound,  which  provides 
greater  stickiness  between  wet  snow  and  the  tire.  We 
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suspect  the  tire  compound  rather  than  the  tread  design 
because  of  the  good  performance  of  the  buffed  (treadless) 
Michelin  in  dry  snow  conditions  (Fig.  lb). 

Military  vehicles 

The  data  obtained  from  the  military  vehicles  are 
given  in  Table  5  and  plotted  on  Figure  2.  Here  again  the 


vehicles  were  tested  on  several  different  types  of  snow 
conditions  and  eq  2  (undisturbed  snow  underlain  by 
packed  snow  or  soil)  is  drawn  on  the  plot  for  reference. 
The  range  of  traction  values  obtained  with  the  CIV  on  ice 
(Blaisdell  and  Harrison  1981 )  with  various  tires  is  also 
shown  in  the  figure.  The  following  observations  can  be 
made  from  Figure  2:  1)  there  is  a  slight  reduction  in 


h.  Buffed  Michelin  XCH4. 

Figure  I .  Traction  data  for  the  CIV  equipped  with  three  different  tire  types. 
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traction  on  hard-packed  snow  (approximately  -2 1  %)  as 
compared  to  undisturbed  snow,  2)  the  traction  obtained 
when  undisturbed  snow  overlays  an  ice  cover  is  signifi- 
cantly  reduced  (approximately  -51%),  and  is  slightly 
higher  than  that  observed  on  ice  alone  and  3)  the  values 
obtained  under  natural  (undisturbed)  snow  conditions 
with  the  HMMWV  and  Bradley  agree  well  with  eq  2 


(i.e.,  the  data  obtained  in  a  prior  field  season).  Least 
squares  regression  analysis  yielded  these  two  equations 
for  the  above  conditions 

Hard-packed  snow:  Tg  =  0.321  N0-91  (7) 

Snow  over  ice:  Tg  =  0.127  /V106.  (8) 


c.  Michelin  XCH4. 
Figure  1  (cont'd). 


Normal  Stress  (kPa) 


Figure  2.  Traction  data  for  the  military  vehicles. 
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Table  5.  Military  vehicles’  traction  data. 


Gross  traction  Normal 

Dale  Coeff.  Stress  stress  Density  S/iotv 

(1989)  ( kPa )  (kPa)  ( kglm 3)  description 


a.  HMMWV 
Weight:  33.450  N 
Contact  area:  0.074  m3 


17  Jan 

0.335 

37.9 

113.0 

180 

undisturbed 

0.294 

33.2 

II  3.0 

560 

hard-packed 

0.222 

25.1 

113.0 

120 

undisturbed  over  ice 

19  Jan 

0.347 

39.2 

113.0 

160 

undisturbed 

20  Jan  0.411  46.4 

b.  HEMTT 

Weight:  268.560  N 

Contact  area  (at  241/276  kPa):  0.149  m2 
Contact  area  (at  138/207  kPa):  0.171  m2 

113.0 

160 

undisturbed 

14  Jan 

0.187 

36.7 

196.3 

250 

undisturbed  over  ice 

0.158 

35.6 

225.3 

250 

undisturbed  over  ice 

18  Jan 

0.282 

63.5 

225.3 

560 

hard-packed 

c.  SUSV 

Weight:  61,340  N 
Contact  area:  1.18  m2 

0.149 

33.6 

225.3 

250 

undisturbed  over  ice 

17  Jan 

0.138 

1.8 

13.0 

140 

undisturbed  over  ice 

0.342 

4.4 

13.0 

560 

hard-packed 

19  Jan 

d.  Bradley 

Weight:  223.299  N 
Contact  area:  2.09  m2 

0.346 

4.5 

13.0 

120 

undisturbed 

14  Jan 

0.182 

9.7 

53.4 

560 

hard-packed 

17  Jan 

e.  M60A1 

Weight:  444,820  N 
Area:  2.8  m2 

0.238 

12.7 

53.4 

230 

undisturbed 

17  Jan 

0.332 

26.4 

79.4 

560 

hard-packed 

At  this  point  one  might  question  the  difference  in 
results  between  the  military  vehicles  and  the  CIV  on 
hard-packed  snow.  The  CIV  (with  the  control  Micheiin 
tires)  obtained  only  slightly  less  traction  than  predicted 
by  eq  2,  and  compares  well  to  data  from  the  other 
vehicles,  as  seen  in  Figure  3.  The  buffed  Micheiin  tires 
generated  more  traction  than  predicted  by  eq  2;  this 
seems  to  indicate  the  effect  of  increased  contact  area  (no 
voids  on  the  tire’s  surface).  The  low  traction  values 
obtained  with  the  NDCC  tire  (as  discussed  above)  are 
seen  to  fall  well  below  the  other  tires. 

Resistance 

The  amount  a  vehicle  sinks  in  the  snow  (sinkage)  or 
the  depth  of  the  rut  left  by  vehicle  passage  greatly  affects 


the  amount  of  motion  resistance  exerted  by  the  snow. 
The  shallow  snow  model  calculates  resistance  based  on 
the  amount  of  sinkage.  The  maximum  sinkage  predicted 
by  SSM2.0  is  determined  from  eq  3.  SSM2.0  uses  Table 
1  to  obtain  the  theoretical  final  density  (pf)  for  the  snow, 
using  nominal  vehicle  contact  pressure  to  enterthe  table. 
We  know  of  no  systematic  study  of  actual  rut  density  for 
wheeled  vehicles,  presumably  because  of  the  difficulty 
in  measuring  densities  of  this  generally  small  layer  of 
dense  snow.  Thus,  it  is  customary  toestimate  or  calculate 
a  compacted  density. 

Snow  properties,  rut  depth  and  densities  were  mea¬ 
sured  forcomparison  with  those  calculated  in  the  SSM2.0. 
Table  6  presents  the  measured  snow  and  rut  characteris¬ 
tics  and  the  sinkage  values  obtained  using  eq  3  where 


8 


Figure  3.  Vehicle  traction  on  hard-packed  snow. 


v 


Figure  4.  Chalk  dust  displacement  around  a  tire  rut. 


final  density  is  obtained  from  Table  1.  In  general,  the 
measured  rut  density  values  do  not  agree  very  well  with 
those  used  in  SSM2.0;  however,  the  values  obtained  for 
sinkage  are  very  close. 

Equation  3  was  derived  by  assuming  that  the  width  of 
the  deformed  snow  is  a  constant  that  is  equal  to  the  width 
of  the  tire  or  track  and  that  the  deformed  snow  is  at  a 
uniform  density.  In  an  effort  to  observe  how  the  snow 
deforms,  we  marked  the  snow  by  filling  1 2.5-mm  (0.5-in.) 
diameter  vertical  holes  punched  in  the  snow,  perpen¬ 


dicular  to  the  direction  of  wheel  travel,  with  chalk  dust. 
A  vehicle  was  then  driven  into  the  marked  area.  Figure 
4  shows  the  results  of  such  an  experiment  and  indicates 
that  there  is  significant  lateral  snow  displacement.  To 
account  for  this  lateral  deformation,  the  deformed  region 
under  the  tire  rut  (Fig.  4)  can  be  assumed  to  have  an 
average  width  h.x  while  the  tire  (track)  has  a  width  of  h. 
Using  the  same  volumetric  procedure  used  to  develop  eq 
3,  we  obtained  the  following 
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Table  6.  Vehicle  sinkage  data. 


Sinkage 


Date 

(1989) 

Snow  Density 
Initial  Final* 

Depth 

(cm) 

Measured 

(cm) 

SSM2.0 

(cm) 

a.  CIV  Tire  A  (179.3-kPa  inflation  pressure) 

13  Jan 

200 

_ 

4.5 

2.0 

2.9 

200 

— 

18.0 

16.5 

11.5 

18  Jan 

90 

— 

7.5 

6.0 

6.3 

19  Jan 

110 

— 

12.0 

10.0 

9.6 

20  Jan 

120 

— 

20.0 

17.0 

15.6 

2  Mar 

150 

— 

20.0 

16.5 

— 

b.  CIV  Tire  D  (248.2-kPa  inflation  pressure) 

14  Jan 

240 

_ 

3.5 

2.5 

2.0 

240 

— 

16.5 

13.0 

9.3 

19  Jan 

110 

— 

10.0 

8.0 

8.0 

110 

— 

13.0 

10.5 

10.4 

125 

— 

8.0 

6.5 

6.2 

20  Jan 

120 

— 

19.0 

14.5 

14.9 

c.  CIV  Tire  C  (413.7-kPa  inflation  pressure) 

19  Jan 

110 

_ 

12.0 

11.0 

10.0 

19  Jan 

95 

— 

6.0 

5.0 

5.1 

20  Jan 

120 

— 

17.0 

16.0 

13.9 

d.  HMMWV  (151.7-kPa  inflation  pressure) 

13  Jan 

200 

380F 

15.0 

8.0 

9.1 

200 

480R 

15.0 

8.0 

9.5 

200 

470F 

19.5 

11.9 

11.8 

200 

500 R 

19.5 

13.8 

12.4 

200 

535F 

23.5 

19.1 

14.3 

200 

— R 

23.5 

21.0 

15.0 

200 

455F 

19.5 

14.5 

11.8 

200 

470R 

17.5 

13.1 

12.4 

200 

445F 

18.0 

14.5 

11.0 

200 

475R 

18.0 

13.0 

11.5 

19  Jan 

120 

320F 

13.0 

9.0 

9.7 

120 

410R 

12.0 

8.5 

10.2 

20  Jan 

120 

— F 

19.5 

14.0 

14.5 

120 

— R 

18.0 

14.0 

15.2 

23  Jan 

190 

575R 

10.0 

9.0 

6.5 

e.  HEMTT  (275.8-kPa  inflation  pressure) 

14  Jan 

250 

490F 

20.0 

12.4 

10.2 

250 

440 

21.0 

16.9 

10.2 

250 

520 

20.0 

16.2 

10.9 

250 

460R 

20.0 

15.6 

10.9 

250 

490F 

16.5 

12.4 

8.4 

250 

490 

18.0 

13.2 

8.4 

250 

510 

20.5 

15.4 

9.0 

250 

510R 

20.0 

15.6 

9.0 

245 

— F 

20.0 

14.3 

10.4 

245 

460R 

17.0 

12.6 

11.1 

20  Jan 

120 

— F 

13.0 

10.0 

9.5 

120 

— R 

13.0 

— 

10.2 

10 


Table  6  (cont’d). 


_ Sinkage _ 

Dale  Snow  Density  Depth  Measured  SSM2.0 

(1989)  Initial  Final*  (cm)  (cm)  (cm) 


f.  HEMTT  (206.8-kPa  inflation  pressure) 


14  Jan 

— 

— R 

6.0 

1.0 

— 

— 

— 

3.5 

3.2 

— 

245 

510F 

17.0 

12.6 

7.7 

245 

525 

18.0 

15.1 

7.7 

245 

430 

17.5 

14.6 

9.4 

245 

445R 

16.0 

13.1 

9.4 

245 

445R 

16.0 

12.2 

8.9 

19  Jan 

— 

— F 

3.5 

0.3 

— 

95 

— R 

4.0 

3.0 

3.3 

120 

300F 

12.0 

10.5 

7.6 

120 

315R 

11.0 

9.0 

9.4 

SUSV  (13.2-kPa  contact  pressure) 

13  Jan 

240 

415F 

14.0 

_ 

7.3 

240 

450R 

15.5 

11.4 

7.3 

250 

460F 

18.5 

10.2 

9.3 

250 

450R 

21.0 

10.8 

9.3 

250 

440F 

17.5 

8.9 

8.8 

250 

430R 

17.0 

9.5 

8.8 

20  Jan 

120 

230F 

12.5 

12.5 

9.5 

120 

245R 

12.5 

12.5 

9.5 

23  Jan 

410 

430R 

64.0 

19.5 

11.5 

265 

455R 

54.0 

27.5 

25.4 

24  Jan 

310 

450R 

58.4 

22.9 

22.2 

450 

480R 

35.5 

5.0 

3.6 

Bradley  (53.09-kPa  contact  pressure) 

14  Jan 

_ 

_ 

3.0 

2.0 

_ 

20  Jan 

100 

355 

18.5 

16.0 

14.8 

23  Jan 

190 

— 

8.0 

7.0 

4.9 

265 

540 

51.0 

37.0 

24.0 

i.  M60A1  (79.4-kPa  contact  pressure) 

16  Jan  340  520  56.0  32.0  17.9 


*F — refers  to  the  front  axle,  R — to  the  rear  axle,  values  between  F  and  R  are  the 
second  and  third  axles  respectively. 


Since  (ft/f»a)  is  less  than  1 ,  the  effect  of  this  term  is  to 
increase  the  sinkage.  This  explains  why,  even  with  the 
estimates  of  final  theoretical  density  (pf)  being  too  high, 
eq  3  still  yields  good  predictions  of  sinkage  in  SSM2.0. 

The  predicted  sinkage  values  in  Table  6  were  deter¬ 
mined  following  the  SSM2.0  algorithm  (eq  3)  with  the 
caveat  that  the  snow  depth  used  for  the  front  tires  was 
also  used  for  all  following  tires,  even  though  the  mea¬ 
sured  snow  depths  (at  the  following  wheels)  may  have 
been  different.  The  use  of  this  algorithm  produces  fairly 
good  results. 


The  measured  sinkage  data  from  different  vehicles 
are  depicted  in  Figure  5.  The  y-axis  variable  (p x/j)  and 
the  A  -axis  variable  (h-z)  were  obtained  by  rearranging  eq 
3  so  that  the  slope  on  this  plot  would  represent  the  final 
density  (pf).  The  data  fall  quite  neatly  on  a  straight  line, 
which  has  a  slope  of  563  kg/m3.  Predictions  of  sinkage 
made  using  the  least-squares  regression  equation 
(y=563.0x+7.1 1)  are  only  slightly  different  (and  no 
better  on  the  average)  from  those  obtained  using  the 
algorithm  of  SSM2.0. 

Motion  resistance  data  for  the  CIV  are  given  in  Table 
7.  Motion  resistance  tests  were  conducted  using  three 
different  width  tires  to  determine  the  effect  of  tire  width 
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on  resistance.  Comparing  the  data  for  tests  conducted 
sequentially  in  the  same  test  area  with  similar  snow 
conditions  revealed  no  clear  trend  when  tire  width  was 
the  only  changing  parameter.  We  wanted  a  parameter 
that  could  be  used  tocombine  vehicle  characteristics  and 
varying  snow  conditions.  After  several  iterations,  the 
parameter  (p 0xbxa)  was  tried,  where  a  is  the  length  of 
the  tire  or  track  in  contact  with  undisturbed  snow  (Fig. 
A 1 ),  and  is  determined  from  the  calculated  sinkage  and 
tire  or  track  geometry.  The  CIV  resistance  attributable  to 
snow  deformation  reported  here  and  that  from  Blaisdell 
et  al.  (1990)  are  plotted  in  Figure  6a.  The  scattered  data 
in  Figure  6a  tend  toward  0  as  the  snow  parameter 
approaches  0.  A  0  intercept  is  expected  since  no  snow 
yields  no  external  resistance.  An  equation  of  the  form 
y=dx?  was  fit  to  the  data  in  Figure  6b. 

Table  8  contains  the  resistance  data  obtained  from  the 
military  vehicles  during  both  field  seasons  (see  Table  3); 
also  included  are  calculations  of  the  parameter  (pQxbxa) 
as  defined  above.  The  resistance  value  presented  is  per 
tire  (track)  and  thus  all  vehicle  data  can  be  shown  on  the 
same  graph.  The  SUSV  was  considered  a  two-track 
vehicle  (as  opposed  to  four),  since  it  appears  that  the 
trailing  track  does  not  act  like  a  trailing  tire  on  a  wheeled 
vehicle  owing  to  its  extremely  low  ground  pressure. 

The  division  of  Rs  by  the  number  of  wheels  or  tracks 
is  a  departure  from  traditional  thoughts  on  vehicle  resis¬ 
tance  in  snow.  In  the  past,  and  in  SSM2.0,  it  was  believed 


that  trailing  tires  with  contact  pressures  lower  than  the 
preceding  tires  would  have  little  or  no  external  resis¬ 
tance.  By  dividing  by  the  number  of  wheels,  we  are 
assuming  here  that  following  tires  displace  as  much 
snow  (on  a  mass  basis)  as  preceding  tires.  This  theory 
was  not  explicitly  examined  during  these  tests,  and  the 
sinkage  measurements  discussed  above  do  not  apply 
since  sinkage  was  measured  behind  each  axle  after  the 
vehicle  came  to  a  stop.  To  test  this  idea,  the  sinkage 
should  be  measured  at  the  same  location  in  the  snow  as 
each  axle  passes  by. 

The  resistance  per  wheel  (per  track)  are  plotted  in 
Figure  7  and  are  combined  in  Figure  8.  The  results  of 
least-squares  regression  analysis  based  on  an  equation  of 
the  form  y=dxe  are  shown  on  each  plot.  The  equation 
obtained  from  using  all  the  data 

Rs  =  68.083  (p0ow)°-9135 

(correlation  coefficient  r2  =  0.39)  ( 1 0) 

is  used  to  estimate  resistance  for  the  vehicles  and  is 
compared  with  both  the  measured  values  and  the  results 
from  SSM2.0  in  Tables  7  and  8. 

Table  9  summarizes  the  percent  errors  shown  in 
Tables  7  and  8.  From  this  table  it  can  be  seen  that  eqlO 
unfortunately  does  not  offer  any  improvement  over 
SSM2.0. 
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*  Percent  difference  from  the  measured  net  resistance. 


Table  8.  Motion  resistance  data  for  military  vehicles. 
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Average:  92  66 

Standard  deviation:  3 1 1  296 


a.  Tracked  vehicles  (cont’d), 
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Table  8  (cont’d).  Motion  resistance  data  for  military  vehicles. 
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Average: 

Standard  deviation: 


b.  Wheeled  vehicles  (cont’d). 
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Resistance  Per  Wheel  (N)  __  Resistance  Per  Wheel  (N) 


h.  Log-log  scale. 


Figure  6.  CIV  resistance  data. 


N,  Net  Resistance  Resistance  Per  Wheel  (N) 


c.  Bradley. 


Figure  7  (cont’d).  Resistance  data  for  military  vehicles. 
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Resistance  Per  Wheel  (N)  _  Resistance  Per  Wheel  (N) 


e.  HEMTT. 


f.  LAV. 

Figure  7  (cont’d). 
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Density  •  Arc  Length  •  Width  (kg/m) 
Figure  8.  Resistance  data  for  all  vehicles. 


Table  9.  Summary  of  percent  differences  between 
resistance  calculations. 


SSM2.0  Eg  10 


Average 

Standard 

deviation 

Average 

Standard 

deviation 

susv 

-52 

-110 

102 

909 

CIV 

36 

156 

60 

98 

Ml  13 

66 

296 

92 

312 

Bradley 

-81 

20 

-74 

32 

HMMWV 

22 

68 

51 

80 

LAV 

-52 

36 

-3 

56 

HEMTT 

-41 

59 

-19 

63 

Avg. 

-15 

75 

30 

221 

Min. 

-81 

-110 

-3 

32 

Max. 

66 

296 

102 

909 

Average  of 

absolute  values:  50  57 


CONCLUSIONS  AND 
RECOMMENDATIONS 

In  this  study  we  have  developed  some  empirical 
equations  that  further  define  gross  vehicle  traction,  spe¬ 
cifically  for  conditions  of  hard-packed  snow  and  for 
snow  over  ice. 

The  sinkage  equation  used  in  the  SSM2.0  was  shown 
to  produce  good  results  and  no  improvement  was  ob¬ 


tained  by  modifying  the  equation  based  on  field  observa¬ 
tions  of  snow  compaction. 

The  data  base  of  resistance  measurements  for  ve¬ 
hicles  in  snow  has  been  increased  and  a  parameter 
developed  that  appears  to  describe  the  wheel-or  track- 
snow  interaction  well.  Vehicle-specific  empirical  analy¬ 
sis  produced  promising  results;  however,  when  all  the 
vehicle  data  are  combined,  an  adequate  predictive  rela¬ 
tionship  is  still  not  obtained.  Vehicle  motion  resistance 
using  SSM2.0  is  only  predicted  to  within  50%  on  aver¬ 
age. 

The  primary  problem  with  resistance  analysis  seems 
to  be  with  trailing  tires  and  how  their  resistance  should 
be  handled.  The  assumption  used  here  that  trailing  tires 
produce  as  much  resistance  as  leading  tires  needs  to  be 
confirmed  experimentally.  Future  shallow  snow  mobil¬ 
ity  work  should  concentrate  in  this  area. 
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APPENDIX  A:  SHALLOW  SNOW  MOBILITY  MODEL  (SSM2.0) 

Description 

The  two  principal  quantities  governing  mobility  are  gross  traction  Tg  and  external  motion  resis¬ 
tance  Rs.  This  model  produces  estimates  of  both  Tg  and  Rs  and  calculates  net  traction  Tn  from  their 
difference 

Tn  =  Ts-Rs.  (Al) 

If  the  net  traction  is  greater  than  0  the  vehicle  is  mobile,  otherwise  it  is  immobilized. 

Motion  resistance  ( Rs )  is  the  resistance  generated  by  terrain  deformation,  in  our  case  by  snow 
compaction.  Compaction  is  partially  the  result  of  vertical  forces  (vehicle  weight)  applied  to  the  snow 
surface  by  the  tire;  however,  it  takes  place  along  a  curved  path  and,  therefore,  horizontal  forces  are 
also  applied.  When  compaction  occurs,  it  can  be  witnessed  by  the  presence  of  a  rut  in  the  snow 
following  vehicle  passage. 

Motion  resistance  Rs  is  a  function  of  many  parameters.  A  partial  list  includes  the  load,  contact 
pressure,  snow  strength  and  depth,  and  width  of  the  tire  or  track.  During  the  past  30  years,  several 
resistance  models  have  been  proposed  in  the  literature  for  deformable  materials.  With  the  goal  of 
keeping  SSM2.0  simple  (i.e.,  a  model  that  has  a  short  list  of  input  parameters  that  are  easily  obtained) 
and  allowing  it  to  address  a  broad  range  of  vehicle  and  snow  conditions,  these  resistance  expressions 
were  scrutinized.  The  vehicle  data  required  by  each  of  these  expressions  are  similar  from  model  to 
model,  and  are  readily  accessible.  The  snow  data  required  to  process  any  of  these  expressions, 
however,  vary  considerably.  The  only  model  that  requires  snow  data  that  can  be  quickly  and  routinely 
obtained  in  the  field  is  that  of  Liston  (1974). 

Liston  assumes  that  a  hyperbolic  relationship  exists  between  compacting  pressure  and  volume. 
Applying  energetics,  he  then  integrates  between  the  initial  and  final  volumes  of  snow  to  obtain  the 
work  of  snow  compaction.  If  no  volume  change  occurs  (no  sinkage  z),  no  work  is  done.  Finally,  the 
work  of  compaction  is  equated  to  external  motion  resistance  Rs  times  the  horizontal  distance  traveled. 

If  we  assume  that  lateral  flow  of  the  snow  during  compaction  is  insignificant  (i.e.,  compaction  is 
confined  to  the  width  of  the  tire  or  track),  volume  change  in  the  snow  can  be  expressed  in  terms  of 
the  sinkage  z.  Further,  if  the  total  mass  of  the  snow  does  not  change  during  compaction,  then  initial 
and  final  volumes  of  snow  can  be  related  to  the  initial  and  final  densities  of  snow.  We  can  then  write 

/ 

Rs  =  pbhp0  ( [  lApf-  po)]  In  (pf/p0)  -  ( 1/Pf) > .  p0  <  pf  (z  >  0)  (A2) 

Rs  =  0  po  =  pf(z  =  0) 

where  p  =  tire  inflation  pressure 

b  =  maximum  tire  or  track  width 
h  =  snow  depth 

p  =  snow  density  after  passage  of  a  given  tire 
p0  =  initial  snow  density  (prior  to  tire  or  track  passage) 
pf  =  maximum  (final)  snow  density  (after  vehicle  passage) 
z  =  sinkage. 

Equation  A2  is  used  in  SSM2.0  for  the  calculation  of  motion  resistance  attributable  to  snow 
compaction. 

Driving  traction  is  also  a  sum  of  the  interaction  of  many  snow  and  vehicle  parameters.  Those  that 
were  mentioned  above  for  resistance  still  apply,  along  with  more  detailed  features  of  the  tire  or  track 
(e.g.,  tread  pattern,  tire  or  track  “rubber”  compound,  grouser  spacing  and  height,  grouser  or  tread 
geometry).  The  number  of  traction  models  proposed  in  the  literature  is  fewer  than  is  the  case  for 
motion  resistance.  These  models  seem  to  fall  into  two  categories,  either  they  are  very  simplistic, 
lumping  many  parameters  together  into  a  few  constants,  or  they  are  exceedingly  detailed. 
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The  SSM2.0  uses  the  Mohr-Coulomb  failure  criteria  relationship  obtained  by  Blaisdell  et  al. 
(1990)  to  determine  gross  traction 

Tg  =  0.851  N0-823  (A3) 

where  N  is  the  normal  stress  under  a  tire  (track)  in  kilopascals.  This  equation  is  based  on  data  from 
a  wide  range  of  vehicles,  but  was  limited  to  initially  undisturbed  snow  conditions  and  for  tracked 
vehicles  or  wheeled  vehicles  equipped  with  tires. 

The  mobility  expressions  given  by  eq  A1-A3  should  be  thought  of  for  a  single  tire  or  traction 
element.  For  a  mobility  model  to  be  flexible,  neither  the  specific  nor  general  configuration  of  a  vehicle 
should  be  limited  by  the  model.  Being  per  tire  or  per  track,  the  relationships  for  traction  and  resistance 
given  here  are  used  on  the  vehicle’s  tires  or  tracks  one  at  a  time  in  SSM2.0  until  all  of  the  traction 
elements  have  been  accounted  for.  In  this  way,  tires  or  tracks  with  different  loads,  inflation  pressures, 
sizes,  configurations  (dual  or  single,  driven  or  free-wheeling),  degree  of  tracking  and  position  on  the 
vehicle  are  all  accommodated  for  with  one  set  of  equations.  Placed  in  a  loop  in  SSM2.0,  these 
expressions  are  used  for  each  station,  and  a  running  sum  for  the  whole  vehicle  is  accumulated  to 
produce  a  measure  of  the  net  performance  of  the  vehicle.  A  station  is  defined  here  as  a  transverse 
section  of  the  vehicle  including  a  single  axle  (single  or  dual  tires)  or  track  loop  (i.e.,  both  sides  of  the 
vehicle  are  assumed  to  be  similar). 

To  apply  the  traction  and  resistance  equations  above  to  even  a  single  tire  or  track,  it  is  necessary 
that  we  calculate  or  measure  several  parameters.  The  determination  of  Tg  requires  only  the  hard 


h - 1 - H 
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contact  area  of  the  tire  or  track  and  the  vehicle 
weight.  To  calculate  Rs,  initial  and  final  densities 
associated  with  the  passage  of  a  tire  or  track  are 
necessary.  In  the  following  paragraphs,  we  first 
describe  how  these  densities  are  determined  and 
then  proceed  to  show  how  the  mobility  equations  are 
adapted  for  application  to  a  whole  vehicle  in  SSM2.0. 

At  each  station,  the  vehicle  parameters  p,  W  (for 
wheeled  vehicles  only ),  l  (for  tracked  vehicles  only), 
b  and  r  of  the  tire  or  track  are  required  (see  Fig  A1 ). 

For  the  snow,  it  is  required  that  we  know  the 
sinkage  of  each  wheel  (track).  To  determine  the 
sinkage,  it  is  assumed  that  compaction  (within  the 
realm  of  loads  that  are  most  common  for  vehicles) 
only  occurs  in  the  vertical  dimension.  First,  we 
assume  that  the  maximum  sinkage  (z^)  that  occurs 
as  the  result  of  vehicle  passage  can  be  calculated 
from 


=  K 


(A4) 


where  h0  =  depth  of  the  undisturbed  snow 

pr  =  maximum  or  final  density  (Yong  and 
Fukue  1977)  in  any  rut  following  ve¬ 
hicle  passage. 

It  is  reasonable  to  believe  that  zmax  and  thus  pf  will 
occur  under  the  tire  or  track  applying  the  largest 
ground  pressure.  pf  is  thus  determined  in  the  SSM2.0 
based  on  the  highest  pressure  (pmax)  exerted  by  any 
tire  or  track  on  the  vehicle.  The  fourmajorcategories 
are  based  on  experience  and  are  defined  in  Table  1 . 


Figure  A3.  Flow  chart  for  shallow  snow  mobility  model. 
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Intermediate  values  of  sinkage  z  for  tires  or  tracks  with  contact  pressures  (p)  less  than  pmax  were 
determined  by  applying  the  ratio  of  (p/pmax)  to  zmax.  Since  the  pressure-sinkage  relationship  is  obvi¬ 
ously  not  linear  for  compressible  snow  (less  than  500  kg/m3)  (Abele  1970),  we  assume  that  a  power 
function  relates  ground  pressure  to  sinkage.  SSM2.0  calculates  sinkage  (referenced  to  the  original 
snow  depth  h0)  for  a  given  station  from 

z  =  ZmaxO’/Pmax)0-5  <A5> 

(see  Mellor  1964,  Fig.  III-34). 

The  sinkage  z,  at  a  given  station  i  on  the  vehicle,  can  then  be  calculated  from 


1-1 

2i  ~  2max  {p  /Pmax)  ^  _  Z  2j 

J- 1 

Pi  >  Pj_,'  Pj_2' 

Po 

(A6) 

Zj=0 

Pi  -  =  Pi-1'  Pi-2' 

•  Po- 

To  calculate  motion  resistance  Rs  we  need  to  know  the  intermediate  values  of  snow  density  as 
compaction  progresses  from  initial  density  p0  to  final  density  pf.  We  have  already  stated  that  pf  is 
associated  with  maximum  sinkage  (z  =  zmax)  and  p0  corresponds  to  a  sinkage  of  z  =  0.  Recalling  eq 
A4,  and  the  assumptions  that  it  is  based  on,  we  can  find  the  density  beneath  a  particular  station  /  from 


Pi  = 


Pi-i 

1-Mi-l) 


(A7) 


where  pi—1  and  hy_\  are  the  snow  density  and  depth  prior  to  the  passage  of  the  tire  or  track  at  station 
i,  and  z;  (eq  A7)  is  the  sinkage  produced  at  the  current  station  (Fig.  A2). 

Lastly,  we  recognize  that  not  all  of  the  tires  or  tracks  on  a  particular  vehicle  may  be  traveling  in 
undisturbed  snow.  Some  tires  or  tracks  may  follow  exactly  in  the  path  of  a  preceding  element,  or  may 
operate  in  undisturbed  snow  (e.g.,  a  narrow  or  wide  trailer  behind  a  vehicle)  or  may  encounter  both 
compacted  and  uncompacted  snow  (e.g.,  dual  tires  following  a  single  tire).  We  need  then  to  account 
for  the  possibility  of  tires  (tracks)  having  some  percentage  (a)  of  their  width  compacting  new  snow 
while  the  remainder  is  traveling  in  a  previously  created  rut.  The  equations  for/?s  and  Tg  are  therefore 
modified  to  become 


Rs.  =  a 


+  (!-«) 


Pi  b,  h0  p0  ( — 1 —  In  -  -U 

\Pi-Po  P°  P>  I. 

pbi  /ij_i  pj_i  ( - 1 - In  -£L  -  -L] 

Vpi -  Pi-1  Pi-1  Pi/ 


for  pi  >  Pm 


(A8) 


R  s;  =  a 


P,b,h0p0[—l —  lnfi-  -1 

Pi-Po  P°  Pi 


for  pi  =  pj_, 


RSl  =  0 


for  pj  <  pi_, 


Tg  =  0.851  /V0-823.  (A9) 

Equations  A8  and  A9  provide  the  essence  of  the  SS  M2.0.  These  equations  are  executed  for  each  station 
of  the  vehicle  and  a  running  sum  for  traction  and  resistance  accumulated.  The  ultimate  ability  of  the 
vehicle  to  move  (total  net  traction  Tn)  is  then  determined  from 
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(A10) 


Tu  =  irgi-/?si 

1=1 


where  n  is  the  total  number  of  stations  on  the  vehicle.  If  Tn  is  positive,  the  vehicle  is  mobile  and  has 
the  capacity  to  accelerate,  climb  hills  or  pull  a  payload  in  proportion  with  the  magnitude  of  Tn.  A  value 
of  0  indicates  impending  immobilization,  and  a  negative  value  of  7"n  predicts  a  definite  no-go  situation. 
A  copy  of  the  SSM2.0,  in  HP  Basic  computer  code  follows,  along  with  the  output  from  sample  runs, 
and  a  flow  chart  is  provided  in  Figure  A3.  Vehicle  data  for  the  SSM2.0  is  in  Table  A 1 . 


Table  Al.  Vehicle  data  for  SSM2.0. 


CIV,  179  kPa  WHEELED 
GVW  =  24696.4  N 

MAXIMUM  GROUND  (INFLATION)  PRESSURE  179.30  kPa 
AVERAGE  HARD  SURFACE  CONTACT  AREA  0.028  m* 


WEIGHT  FER  TIRE  NEW  SNOW 

STATION  RADIUS  (cm)  WIDTH  (an)  WHEEL  (N)  IRES.  (kPa)  DRIVEN  DUALS  PERCENT 


1  35.6  25.4  6174.1  179.3  Y  N  100 

2  35.6  25.4  6174.1  179.3  N  N  0 


CIV,  103  kPa  WHEELED 
GVW  =  24696.4  N 

MAXIMUM  GROUND  (INFLATION)  PRESSURE  110.30  kPa 
AVERAGE  HARD  SURFACE  CONTACT  AREA  0.0345  Vfi 

WEIGHT  PER  TIRE  NEW  SNOW 

STATION  RADIUS  (cm)  WIDTH  (cm)  WHEEL  (N)  FEES.  (kPa)  DRIVEN  DUALS  PERCENT 


1  34.9  26.7  6174.1  110.3  Y  N  100 

2  34.9  26.7  6174.1  110.3  N  N  0 


HEMTT,  207/276  kPa  WHEELED 
GVW  =  268560.1  N 

MAXIMUM  GROUND  ( INFLATION)  PRESSURE  275.80  kPa 
AVERAGE  HARD  SURFACE  CONTACT  AREA  0.149  ffi* 


STATION 

RADIUS  (cm) 

WIDIH  (cm) 

WEIGHT  PER 
WHEEL  (N) 

TIRE 

IRES.  (kPa)  DRIVEN 

DUALS 

NEW  SNOW 
PERCENT 

1 

61.7 

46.5 

31004.0 

241.3 

Y 

N 

100 

2 

61.7 

46.5 

31459.9 

241.3 

Y 

N 

0 

3 

61.7 

47.5 

35785.8 

275.8 

Y 

N 

0 

4 

61.7 

47.5 

36030.4 

275.8 

Y 

N 

0 
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Table  A1  (cont’d).  Vehicle  data  for  SSM2.0. 

HEfflT,  138/207  kPa  WHEELED 
GVW  =  268560.1  N 

MAXIMUM  GROUND  (INFLATION)  PRESSURE  206.80  kPa 
AVERAGE  HARD  SURFACE  CONTACT  AREA  0.171  m2 


WEKWT  PER  TIRE  NEW  SNOW 

STATION  RADIUS  (an)  WIDTH  (an)  WHEEL  (N)  IRES.  (kPa)  [RIVEN  DUALS  PERCENT 


1 

58.9 

48.9 

31004.0 

137.9 

Y 

N 

100 

2 

58.9 

48.9 

31459.9 

137.9 

Y 

N 

0 

3 

58.9 

48.3 

35785.8 

206.8 

Y 

N 

0 

4 

58.9 

48.3 

36030.4 

206.8 

Y 

N 

0 

HMMWV,  138/152  kPa  WHEELED 
GVW  =  33450.5  N 

MAXIMUM  GROUND  (INFLATION)  PRESSURE  151.70  kPa 
AVERAGE  HARD  SURFACE  OCNEACT  AREA  0.074  m2 

WEIGHT  PER  TIRE  NEW  SNOW 

STATION  RADIUS  (an)  WIDTH  (an)  WHEEL  (N)  PRES.  (kPa)  [RIVEN  DUALS  PERCENT 


1  42.9  32.0  7228.3  137.9  Y  N  100 

2  42.9  33.0  9496.9  151.7  Y  N  0 


LAV, (12.5X20)  207  kPa  WHEELED 
GVW  =  125483.7  N 

MAXIMUM  GROUND  ( INFLATION)  PRESSURE  206.80  kPa 
AVERAGE  HARD  SURFACE  CONTACT  AREA  0.100  ffi2 


WEIGHT  PER  TIRE  NEW  SNOW 

STATION  RADIUS  (an)  WIDTH  (an)  WHEEL  (N)  PRES.  (kPa)  [RIVEN  DUALS  PERCENT 


1 

44.5 

35.3 

16102.5 

206.8 

Y 

N 

100 

2 

44.5 

35.3 

17281.3 

206.8 

Y 

N 

0 

3 

44.5 

35.3 

14501.1 

206.8 

Y 

N 

0 

4 

44.5 

35.3 

14857.0 

206.8 

Y 

N 

0 

LAV, 

(12.5X20)  103  kPa 

WHEELED 

GVW 

=  125483.7  N 

MAXIMUM  GROUND  (INFLATION)  PRESSURE  103.40  kPa 

AVERAGE  HARD  SURFACE  CONTACT  AREA 

0.141  m2 

WEIGHT  PER 

TIRE 

NEW  SNOW 

STATION  RADIUS  (an) 

WIDIH  (on) 

WHEEL  (N)  PRES.  (kPa) 

[RIVEN 

DUALS 

PERCENT 

1 

41.4 

37.8 

16102.5 

103.4 

Y 

N 

100 

2 

41.4 

37.8 

17281.3 

103.4 

Y 

N 

0 

3 

41.4 

37.8 

14501.1 

103.4 

Y 

N 

0 

4 

41.4 

37.8 

14857.0 

103.4 

Y 

N 

0 

LAV,  (11X16)  289  kPa  WHEELED 
GVW  -  119634.3  N 

MAXIMUM  GROUND  (INFLATION)  PRESSURE  289.60  kPa 
AVERAGE  HARD  SURFACE  OCNEACT  AREA  0.058  m2 

WEIGHT  PER  TIRE  NEW  SNOW 

STATION  RADIUS  (an)  WIDTH  (an)  WHEEL  (N)  PRES.  (kPa)  [RIVEN  DUALS  PERCENT 


1 

43.4 

31.4 

15346.3 

289.6 

Y 

N 

100 

2 

43.4 

31.4 

16480.6 

289.6 

Y 

N 

0 

3 

43.4 

31.4 

13822.8 

289.6 

Y 

N 

0 

4 

43.4 

31.4 

14167.5 

289.6 

Y 

N 

0 
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Table  A1  (cont’d), 


LAV, (11X16)  165  kPa  WHEELED 
GVW  =  119634.3  N 

MAXIMUM  GROUND  (INFLATION)  PRESSURE  165.50  kPa 
AVERAGE  HARD  SURFACE  CONTACT  AREA  0.102  in1 2 3 


STATION 

RADIUS  (an) 

WIDIH  (an) 

WEKHT  PER 
WHEEL  (N) 

TIRE 

PRES.  (kPa) 

DRIVEN 

DUALS 

NEW  SNOW 
PERCENT 

1 

41.7 

33.2 

15346.3 

165.5 

Y 

N 

100 

2 

41.7 

33.2 

16480.6 

165.5 

Y 

N 

0 

3 

41.7 

33.2 

13822.8 

165.5 

Y 

N 

0 

4 

41.7 

33.2 

14167.5 

165.5 

Y 

N 

0 

SUSV  TRACKED 
GVW  =  61340.7  N 

MAXIMUM  GROUND  (INFLATION)  PRESSURE  13.20  kPa 
AVERAGE  HARD  SURFACE  CONTACT  AREA  1.180  m2 

WEIGHT  PER  TRACK  NEW  SNOW 

STATION  RADIUS  (can)  WIDTH  (an)  TRACK  (N)  LENGTH  (an)  PERCENT 


1  26.40  60.96  15390.80  190.50  100.00 

2  26.40  60.96  15279.60  190.50  0.00 


Ml  13  TRACKED 
GVW  =  104087.9  N 

MAXIMUM  GROUND  (INFLATION)  FRESSUKE  51.02  kPB 
AVERAGE  HARD  SURFACE  CONTACT  AREA  1.020  m2 

WEIGHT  PER  TRACK  NEW  SNOW 

STATION  RADIUS  (an)  WIDIH  (an)  TRACK  (N)  LENGTH  (an)  PERCENT 


1  36.80  38.10  52043.50  266.70  100.00 


BRADLEY  TRACKED 
GVW  =  223299.6  N 

MAXIMUM  GROUND  (INFLATION)  PRESSURE  53.09  kPa 
AVERAGE  HARD  SURFACE  CONTACT  AREA  2.090  m2 

WEIGHT  PER  TRACK  NEW  SNOW 

STATION  RADIUS  (an)  WIDIH  (an)  TRACK  (N)  LENGTH  (an)  PERCENT 


1  35.60  53.34  111649.75  391.20  100.00 


5-TON,  207  kPa  WHEELED 
GVW  -  105778.2  N 

MAXIMUM  GROUND  (INFLATION)  PRESSURE  206.80  kPa 
AVERAGE  HARD  SURFACE  CONTACT  AREA  0.171  m2 

WEIGHT  PER  TIRE  NEW  SNOW 

STATION  RADIUS  (an)  WIDIH  (an)  WHEEL  (N)  PRES.  (kPa)  DRIVEN  DUALS  PERCENT 


1  55.4  40.8  23842.4  206.8  Y  N  100 

2  55.4  39.2  14523.4  206.8  Y  N  0 

3  55.4  39.2  14523.4  206.8  Y  N  0 
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Computer  code 


10  ! 

20  I  SSM/2.01 

30  I  (trio  tr  sinkage  and  snou  density  function  of  cont .  pr 

40  I  (manual  input  of  all  vehicle  and  snou  data) 

50  I  (Liston's  resist  best-fit  pouer  curve  for  traction) 

60  I 

70  net-0 
80  RR-0 
90  DB-0 
100  go-0 
1 10  pl-0 
120  z-0 
130  sumz-0 

140  DISP  “vehicle?'1 
150  INPUT  UehS 

160  DISP  “wheeled  <w)  or  tracked  (t)?“ 

170  INPUT  type* 

180  DISP  “highest  station  ground  pressure  (kPa)-" 

190  INPUT  pmax 

200  pmax-pmax»0. 1  1  convert  from  kPa  to  N/cm1^ 

210  ! 

220  !  •  •  establish  final  density  based  on  largest  footprint  pressure  *• 

230  sigmaf-0.5 

240  IF  pmax>21  THEN  sigmaf-0.55 
250  IF  pmax>35  THEN  slgmaf-0.6 
260  IF  pmax>70  THEN  sigmaf-0.65 
270  1 

280  DISP  “snou  depth  (cm)  -" 

290  INPUT  h0 

300  DISP  "snow  density  (g/cn“3)  •  “ 

310  INPUT  sigma0 
320  signal -sign*? 

330  sigma2-signa0 
340  h1-h0 

350  1 

360  PRINT  ‘Vehicle:  “lUeh* 

370  PRINT  ‘  * 

380  PRINT  USIN6  formatl  i  “initial  state:  snow  depth  »"ih0i"  cm" 

390  PRINT  USIN6  format2  I  “  snou  density  «"isigma0i"  g/cn*3 

400  ! 

410  zmax-h0»( 1-sigma0/sigmaf )  !  calculate  maximum  sinkage 
420  IF  zmax<0  THEN  zmax-0 
430  ! 

440  RAD  !  compute  in  radians  units  for  trig  functions 
450  I 
460  ! 

470  !  ••  enter  vehicle  data  one  tire  station  at  a  time  •* 

480  DISP  "how  many  wheel  or  track  stations  on  each  side  of  the  vehicle?" 
490  INPUT  N 
500  ! 

510  FOR  1-1  TO  N 
520  DISP  "" 

530  DISP  “station  *jl 

540  DISP  "wheel  radius  or  approx  radius  of  compacting  portion  of  track  cm) 
550  INPUT  rads 

560  OISP  "single  tire  or  track  width  at  this  location  (cm)  -“ 

570  INPUT  wid 

580  DISP  "single  tire  or  track  load  (N>  «" 

590  INPUT  loa 
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600  DISP  “contact  area  (mA2)  *" 

610  INPUT  area 

620  IF  type*-“t“  fl'EN  GOTO  skipl 
630  OISP  "inflation  pressure  (kPa)  =  ” 

640  INPUT  pres 

650  pres-pres*0 . I  1  convert  from  kPa  to  N/cm*2 
660  GOTO  skip2 

670  skipl:  pres-loa/area/ 1 0000  1  convert  from  N/n"2  to  N/cm"2 
680  i 

690  skip2:  OISP  “driven?” 

700  INPUT  rS 

710  IF  type$-"t“  THEN  GOTO  skip3 
720  OISP  “duals  ?* 

730  INPUT  R$ 

740  IF  R$<> “ y "  AND  R$<>”Y"  THEN  GOTO  780 
750  wid-uid»2 
760  loa»loa*2 
770  ! 

780  skip3:  OISP  "percent  of  width  compacting  virgin  snow  (X ) 

790  INPUT  prcnt 
800  I 

810  IF  pres>pl  AND  slgma2<sigmaf  AND  h0>0  THEN  GOTO  870  !  added  sink,  here? 

820  z-0 

830  tempres-0 

840  IF  prcnt>0  AND  sigma2>sigma0  THEN  GOTO  960 
850  resis-0 
860  GOTO  1010 

870  z-<pres/pmax KO.S-zmax-sumz  !  calculate  additional  sinkage  this  station 
880  I 

890  !  «*  set  rut  bottom  values  •• 

900  h2«h1-z 

910  sigma2»sigma1 /< 1 -z/hl ) 

920  ! 

930  !  •*  calculate  resistance  parameter  •• 

940  tempres-1 /( sigma2-sigmal )»L0S< sigmaZ/sigma I  >-l /slgma2  !  in  rut 
950  ! 

960  !  *•  calculate  motion  resistance  at  this  station  •• 

970  resis-prcnt/ 100»< pres»wld«h0»sigma0 >  I  in  virgin  snow 

980  resls»resis*< 1 /< sigmp2-8igma0 )»L06( sigma2/sigma0 >-l /sigma2 >  (virgin  snow 
990  moreres”tempres»< 1 -prcnt/ 1 00 )•< pres«wid#hl ♦slgmal  >  '  in  rut 
1000  resis*resls+moreres 
1010  I 

1020  !  ••  calculate  gross  traction  at  this  station  •• 

1030  trac»0. 851 *< loa/area/ 1000  )“0. 823  I  in  kPa 
1040  trac*trac*area# 1 000  1  in  N 
1050  ! 

1060  IF  r$-"n"  OR  r$="N"  THEN  trac-0  <  no  traction  if  not  driven 
1070  ! 

1080  i  ••  double  to  account  fcr  both  sides  of  the  vehicle  •• 

1090  arrea-2-area 
1100  trac-2#trac 
1110  resis*2»resis 
1 120  ! 

1130  !  ••  print  output  for  this  station  •• 

1140  PRINT  “  " 

1150  PRINT  “station  "«I 

1160  PRINT  USING  format!  i  '  additional  sinkage  «Hizt“  cm" 

1170  PRINT  USING  format!  i  "  total  area-” larreai”  m“2" 

1180  PRINT  USING  format!  i  “  snow  resis  -"tresisi*  N" 

1190  PRINT  USING  formatl  I  "  snow  trac  -“itract“  N“ 
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1200  i 

1210  !  ••  running  summation  for  whole  vehicle  •* 

1220  RR*RR+resis  1  sum  for  whole  vehicle 
1230  0B*0B+trac  !  sum  for  whole  vehicle 

1240  net*trac-resis  1  calculate  net  traction  for  individual  station 
1250  PRINT  USING  format!  i  "  net  snow  traction  for  station  *  "ineti"  N” 

1260  PRINT  USING  format!  i  *  rut  bottom:  depth*" ih2i "  cm” 

1270  PRINT  USING  format2  i  "  density*" isigma2i "  g/cmA3" 

1280  go*go+net  !  sum  net  traction  for  vehicle 
1290  ! 

1300  i  ••  save  last  station  values  for  next  iteration  *• 

1310  IF  pres>pl  THEN  pl=pres 

1320  sigma1*sigma2 

1330  hi *h2 

1340  sumz*sumz+z 

1350  ! 

1360  NEXT  I 
1370  ! 

1 380  !  •»  calculate  mobility  in  English  units  for  output  »• 

1390  eRR*RR/4. 448222 
1400  eOB-OB/4. 448222 
1410  ego«go/4. 448222 
1420  ez*sumz/2.54 
1430  ! 

1440  !  *•  print  out  whole  vehicle  results  •• 

1450  PRINT  *  ' 

1460  PRINT  ’  * 

1470  PRINT  USIN6  format  1  t  “  total  sinkage  for  vehicle*  "jsumzi"  cm  (“jezi*  in 
>“ 

1480  PRINT  USING  farmatl  t  ‘  total  snown  resistance  *  "|RR»“  N  ("ieRRi“  lb>* 
1490  PRINT  USING  format!  i  *  total  snow  traction  *"iDBi“  N  ('leOBi*  1 b  >  “ 

1500  PRINT  "  * 

1510  PRINT  USING  format)  i  ‘  net  traction  for  vehicle  *  "igop*  N  <*iegot*  lb)" 
1520  PRINT  "  " 

1530  PRINT  "  * 

1540  PRINT  "  " 

1550  ! 

1560  format!:  IMAGE  3<K,DDDDD0. DO) 

1570  format2 :  IMAGE  3<K,0Z.0000> 

1580  ! 

1590  END 
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Vehicle:  HMMUV 


initiel  state:  snow  depth  *  5.00  cm 

snow  density  •  0.5S00  g/cmA3 


station  t 

additional  sinkage  *  0.00  cm 

total  area*  0.1480  mA2 
snow  res  is  *  0.00  N 

snow  trac  =  5467.65  N 

net  snow  traction  for  station  *  5467.65  N 

rut  bottom:  depth*  0.00  cm 

density*  0.5500  g/cmA3 


station  2 

additional  sinkage  *  0.00  cm 

total  area*  0.1480  mA2 
snow  res.s  *  0.00  N 

snow  trac  *  6844.84  N 

net  snow  traction  for  station  =  6844.84  N 

rut  bottom:  depth*  0.00  cm 

density*  0.5500  g/cmA3 


total  sinkage  for  vehicle*  0.00  cm  (  0.00  in) 

total  snown  resistance  “  0.00  N  <  0.00  lb) 

total  snow  traction  *  12312.49  N  <  2767.96  lb) 

net  traction  for  vehicle  *  12312.49  N  <  2767.96  lb) 


Vehicle:  HEMTT  (241/276) 

initial  state:  snow  depth  *  25.00  cm 

snow  density  *  0.2200  g/cmA3 


station  1 

additional  sinkage  *  14.03  cm 

total  area*  0.2980  mA2 

snow  resis  *  11515.55  N 
snow  trac  *  20514.00  N 
net  snow  traction  for  station  *  8998.45  N 

rut  bottom:  depth*  10,97  cm 

density*  0.5014  g/cm'3 

station  2 

additional  sinkage  *  0.00  cm 

total  area*  0.2980  mA2 

snow  resis  *  0.00  N 

snow  trac  *  20761.93  N 
net  snow  traction  for  station  ■  20761.93  N 

rut  bottom:  depth*  10.97  cm 

density*  0.5014  g/cmA3 

station  3 

additional  sinkage  ■  .97  cm 

total  area*  0.2980  m*2 

snow  resis  ■  1230.94  N 

snow  trac  ■  23084.34  N 
net  snow  traction  for  station  -  21853.40  N 

rut  bottom:  depth*  10.00  cm 

density*  0.5500  g/cn“3 
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Vehicle:  HEMTT  (24I/27G)  (cont'd). 

initial  state:  snow  depth  -  25.00  cw 

snow  density  *  0.2200  g/cmA3 


station  4 

additional  sinkage  ■  0.00  cm 

total  area*  0.2980  m “ 2 
snow  res  is  •  0.00  N 

snow  trac  -  23214.12  N 
net  snow  traction  for  station  *  23214.12  N 

rut  bottom:  depth"  10.00  cm 

density"  0.5500  g/cm'3 


total  sinkage  for  vehicle-  15.00  cm  (  5.91  in) 

total  snown  resistance  -  12746.49  N  (  2865.52  lb) 

total  snow  traction  *  07574.39  N  (  19687.50  lb) 

net  traction  for  vehicle  *  74827.90  N  (  16821.98  lb) 


Vehicle:  SUSV 

initial  state:  snow  depth  ■  9.00  cm 

snow  density  •  0.1480  g/cmA3 


station  1 

additional  sinkage  -  6.30  cm 

total  area-  2.36  mA2 
snow  resis  -  309.13  N 

snow  trac  -  16626.42  N 
net  snow  traction  for  station  *  16317.29  N 

rut  bottom:  depth"  2.70  cm 

density"  0.4930  g/cmA3 


station  2 

additional  sinkage  -  0.00  cm 

total  area"  2.36  mA2 
snow  resis  -  0.00  N 

snow  trac  «  16527.49  N 
net  snow  traction  for  station  •  16527.49  N 

rut  bottom:  depth-  2.70  cm 

density-  0.4930  g/cmA3 


total  sinkage  for  vehicle-  6.30  cm  (  2.48  in) 

total  snown  resistance  *  309.13  N  (  69.49  lb) 

total  snow  traction  «  33153.90  N  <  7453.29  lb) 

net  traction  for  vehicle  -  32844.78  N  (  7383.80  lb) 
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Figure  B2.  Tire  A  ( Michelin  LT235XCH4; 
radial,  all-season;  179-kPa  inflation  pres¬ 
sure;  0.04 12-m i2  contact  area). 


Figure  B5.  LAV-25  tire  (Michelin  J 2. 5/75  R20  XL;  207-kPa 
inflation  pressure;  0.1 -m2  contact  area). 


Figure  B6.  HEMTT  tire  (Michelin  16.0 
R20;276-lcPa  inflation  pressure;  0.149- 
m2  contact  area). 
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